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Tracking by the GRAIL spacecraft has yielded a model of the gravitational field of the 
Orientale basin at 3—-5-km horizontal resolution. The diameter of the basin excavation 
cavity closely matches that of the Inner Depression. A volume of at least (3.40.2)x10° km® 
of crustal material was removed and redistributed during basin formation; the outer edges 
of the zone of uplifted mantle slope downward and outward by 207-25°. There is no 
preserved evidence of the transient crater that would reveal the basin’s maximum volume, 
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but its diameter may now be calculated from the observed structure to be between the 
diameters of the Inner Depression and Inner Rook ring. The model resolves distinctive 
structures of Orientale’s three rings, including their azimuthal variations, and suggests the 
presence of faults that penetrate the crust. The crustal structure of Orientale provides 
constraints in the third dimension on models for the formation of multi-ring basins. 


Basin-scale impacts disrupted early planetary crusts and imparted substantial kinetic 
energy that was manifest in fracturing, melting, vaporization, seismic shaking, ejection, 
and re-distribution of crustal and possibly mantle material. Impact basins preserve the 
record of these planet-altering events, and the study of these structures elucidates the 
partitioning of energy and its corresponding geological and environmental effects early in 
planetary history. Surface signatures of impact basins on solid planets have been 
extensively documented (/, 2), but their subsurface structure has, to date, been poorly 
characterized. Here we present the highest-resolution orbital gravity field model of any 
major impact basin — the Orientale basin on the Moon — as mapped by the Gravity 
Recovery and Interior Laboratory (GRAIL) mission (3). 

Orientale, located on the western limb of the lunar nearside, is the youngest (~ 3.8- 
Gy old) (/, 4) large (~930-km-diameter) impact basin on the Moon. As a consequence of 
its good state of preservation (/, 5), with relatively few superposed large craters (6), it is 
often considered the standard example of a well-preserved, multi-ring basin in 
comparative studies of large impacts on terrestrial planetary bodies (2, 7). Because of the 
basin’s importance, the GRAIL Extended Mission (XM; see SOM) featured a low- 
altitude mapping campaign during the mission’s Endgame phase (8), in which the dual 
spacecraft orbited the Moon at an average altitude of 6 km and acquired observations less 
than 2 km above the basin’s eastern rings (Figs. $1 and $2). 

To realize the highest-resolution gravity map achievable from the data and to assure 
that small-scale features resolved were robust, we developed two maps that used the same 
data but independent methodologies (see SOM). The first is a global spherical harmonic 
expansion of GRAIL’s Ka-band range-rate (KBRR) tracking data to degree and order 
1200 (spatial block size = 4.5 km). The second is a local model that implemented a short- 
arc analysis (9) of the tracking data and used a gravitational field model to degree and 
order 900 (/0) as the a priori field. Local gravitational anomalies were estimated with 
respect to the spherical harmonic model at the center coordinates of the grid points. The 
final model has a mixed grid resolution 0.1° x 0.1° and 0.1667° x 0.1667°, with 
maximum spatial resolution varying between 3 and 5 km. These independent analyses 
produced gravitational models of Orientale that are essentially indistinguishable (cf. Fig. 
S4). The maps in Fig. 1 resolve the shallow subsurface structure of Orientale at a spatial 
resolution comparable with that of many geological structures at the surface, including 
simple and secondary craters. 

The topography of the Orientale basin from the Lunar Orbiter Laser Altimeter (//) 
and the free-air gravity anomaly field of the region are shown in Figs. 1A and 1B. The 
maps show similar detail at small spatial scales because above degree and order 80 
(spatial block size < 68 km) >98% of the lunar gravity field is attributable to topography 
(12). The high correlation of topography and gravity at short horizontal scales is due to 
the large magnitude of the gravity anomalies arising from topography relative to the 


weaker anomalies arising from density anomalies in the shallow subsurface (/2). 
Removal of the gravitational signature associated with topography clearly resolves 
structure within the lunar crust. 

Both topography and free-air gravity anomaly resolve Orientale’s Inner Depression 
(ID), as well as the Inner Rook ring (IRR), Outer Rook ring (ORR), and Cordillera ring 
(CR). The rings, partially resolved in pre-GRAIL gravitational models with an order of 
magnitude poorer spatial resolution and multiple orders of magnitude lower signal 
sensitivity than those in this study (/3), formed in the process of cavity collapse during 
the modification stage of the impact event, within an hour of the initiation of basin 
formation (/4). The mechanism for ring formation, however, has been controversial (2, 5, 
15-17), in large part because of a lack of understanding of subsurface structure needed to 
provide constraints on impact basin formation models. 

At larger spatial scales, a greater fraction of the gravitational signal is associated with 
mass variations within the Moon. These variations are best revealed in the Bouguer 
gravity anomaly field (Fig. 1C), a representation of the gravitational field after the 
attraction of surface topography has been removed. 

For uniform-density crust and mantle (2550 and 3220 kg m®, respectively), the 
Bouguer gravity anomaly can be downward continued and used to map the crust-mantle 
boundary, to yield crustal thickness (Fig. 1D, Fig. S5). Although the assumption of 
uniform density is an approximation, its application to the regional crustal structure is 
supported by crustal density inferred from GRAIL (/8) as well as from orbital remote 
sensing data (see SOM and Fig. S6). Some models for the crust invoke a mixed 
feldspathic layer that overlies a layer of pure anorthosite (/9), but the density contrast 
between these rock types is small in comparison with that across the crust-mantle 
interface. 

A cross-section of crustal structure along the profile in Fig. 1D is shown in Fig. 2. 
Crust-mantle boundary profiles take into account the effect on the crustal structure of a 
10-km-thick sheet of density 2650 kg m” appropriate for crustal impact melt or crustal 
material formed by differentiation of Orientale’s melt pool (see SOM); in the end- 
member case in which the sheet is instead identical in density to feldspathic crust, the 
minimum crustal thickness would be only 2 km less. Consequently, the presence of a 
high-density sheet does not have a substantial effect on the basin’s crustal structure. Fig. 
1D indicates that the Orientale impact removed and redistributed a minimum of (3.4+0.2) 
x 10° km? of material (see SOM) from the lunar crust. Approximately one-third of the 
excavated material was deposited as ejecta in an annulus between the Cordillera ring and 
a radial distance from the basin center of one basin diameter (20), contributing to 
enhanced crustal thickness in this region (Fig. 1D). The high porosity (~18%) observed in 
the Orientale ejecta blanket (/8, 2/) is similar to the porosity observed in Apollo 14 
samples of the Fra Mauro Formation, which consist of ejecta from the Imbrium basin 
(22). 

The calculation of crustal thickness does not include the presence of the mare units 
that are generally less than 1 km thick and irregularly distributed within the Inner 
Depression and ponded irregularly inside the Outer Rook and Cordillera rings (27) (low- 
reflectance areas in Fig. S7). The gravitational signal of the mare units is partially masked 
by the unknown characteristics of the impact melt sheet, and where these deposits are 
present, the underlying crustal thickness would be slightly greater than in the model. 


Pre-GRAIL gravitational models have shown that the most distinct, well-preserved, 
large-scale gravitational signature associated with large lunar basins is a strong central 
positive free-air anomaly (23-26). These mass excesses or mascons (27) are a 
consequence of the thinning of the crust beneath the basin combined with the super- 
isostatic uplift of a central basin that was mechanically coupled to initially sub-isostatic 
surrounding material (28, 29). Figs. 1 and 2 show that in Orientale, the diameter of 
crustal thinning corresponds closely to the diameter of the Inner Depression. We 
therefore identify the Inner depression as an approximation to the basin excavation 
cavity. 

Of great interest in impact studies is the transient crater — the cavity formed at the 
culmination of excavation, prior to collapse and modification. The transient crater 
represents the time corresponding to the maximum displaced volume, and its diameter 
satisfies standard impact scaling laws (30, 31), which allow many characteristics of an 
impact to be inferred from its size. From previous analyses ranging from geological 
mapping to numerical simulation, each of Orientale’s basin rings have been variously 
interpreted as marking the diameter of the transient crater (2, 5, 15, 16, 32). GRAIL’s 
greatly improved spatial resolution and signal sensitivity compared with all previous 
models of the internal structure of Orientale, in contrast, show no evidence of the 
transient crater, indicating that any remnant is not preserved or at least not recognizable 
in the gravity field. However, the transient crater diameter may now be calculated with 
confidence from scaling laws (33). From the 530-km diameter of crustal thinning (Fig. 
1D) measured from an azimuthally averaged crustal thickness profile (Fig. $5), the 
transient crater diameter is between 320 and 460 km, placing it between the diameters of 
Orientale’s Inner Depression and Inner Rook ring. The transient crater thus does not 
correlate with a specific ring; indeed, hydrocode modeling constrained by this crustal 
structure model (/4) indicates that rings form subsequent to the transient crater, during 
the collapse phase. 

The sharp transition between the basin excavation cavity and the surrounding crust is 
well illustrated in Fig. 2. At the outer edges of the zone of mantle uplift, the crust-mantle 
boundary slopes outward and downward by 20°—25°. The spatial correspondence of this 
plug of uplifted mantle with the Inner Depression is similar to the pattern seen in other 
multi-ring basins (34), but is in contrast to peak-ring basins, where the zone of uplifted 
mantle is limited to within the peak ring (35). 

The model also shows, beyond the basin depression, an annulus of thickened crust 
(Figs. 1D and S8D) and radial structure in gravity gradients (Figs. 1E and S8E) that 
correlate with observed ejecta structures (e.g., secondary crater chains) (20). 

Aspects of Orientale’s asymmetry in surface structure extend to the subsurface, as 
indicated in Figs. 1 and 2. For instance, the basin exhibits an east—west variation in 
regional crustal structure that predated formation of the basin. There are also azimuthal 
variations in crustal thickness, some of which are distinctly associated with basin rings. 
The simplest interpretation of the azimuthally averaged models is that they could 
correspond to displacements associated with normal faults that penetrate the crust. The 
crust-mantle boundary relief in Fig. 2B suggests there could be multiple faults dipping 
inward from the Outer Rook and Cordillera rings. The crustal thickness model also 
suggests the presence of other crustal faults that lack a visible surface expression. 
Although these faults may be listric, 1.e., the dip angle decreases with depth, a dip of 50°, 


indicated by hydrocode simulations (/4), is consistent with prominent changes in crust— 
mantle boundary depth. The simulations also support crustal faulting not associated with 
rings. 

Insight into the distinctive nature of each ring can be gained from scrutiny of Figs. 1, 
2, and S8. The Inner Depression has the most axisymmetric and the largest variation in 
crustal thickness as well as a change in the sign of the gravity gradient. The Inner Rook 
ring has a well-developed topographic expression. Individual peaks within the Inner 
Rook are associated with positive free-air and Bouguer anomalies embedded within an 
annulus of negative free-air and Bouguer anomalies. The Inner Rook also appears 
associated with a near-circumferentially continuous change in the sign of the gravity 
gradient (Fig. 1E) and a flattening in relief along the crust-mantle boundary. 

The Outer Rook ring also displays well-developed topography and associated positive 
free-air anomalies embedded within the same annulus of negative free-air and Bouguer 
anomalies. The most negative Bouguer gravity in the region appears within the Outer 
Rook and may reflect a combination of thickening of the crust by ejecta and extensive 
fracturing in the crustal column. The Outer Rook displays a sign change in the gravity 
gradient and a mild shoaling of the crust-mantle boundary. These characteristics of the 
Outer Rook ring are also consistent with local thinning of the crust associated with 
faulting. 

The topography of the Cordillera ring deviates markedly from axisymmetry; it is less 
developed than the Inner and Outer Rook rings and has little expression in part of the 
basin’s southwestern quadrant. This asymmetric structure may be a consequence of the 
northeast-to-southwest-directed oblique impact that formed the basin (36) or pre-existing 
heterogeneity of crustal or lithospheric structure (2, 37), with a clear west-to-east gradient 
of decreasing crustal thickness still preserved (Fig. 1D). The ring is characterized by 
positive free-air and Bouguer anomalies, a gradient in crustal thickness, and a 
circumferentially discontinuous sign change in the gravity gradient. The variation of 
relief along the crust—mantle boundary strongly suggests fault penetration to the lower 
crust and possibly upper mantle. The gravitational signature could alternatively reflect 
contributions from magmatic intrusions along the irregularly developed ring fault, but 
regional seismic reflection profiles of a portion of the terrestrial Chicxulub impact 
structure, 20%-25% the size of Orientale, show ring faults that extend well into the lower 
crust (38). 

A hydrocode calculation (/4) constrained by the high-resolution observations 
reported here successfully reproduces the major characteristics of Orientale’s crustal 
structure and elucidates the thermal, tectonic, and geological consequences of the largest 
impacts that dominated the earliest evolution of the Moon and other solid planets. 
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FIGURE CAPTIONS AND FIGURES 


Fig. 1. 


Fig. 2. 


(A) Topography, (B) free-air anomaly, (C) Bouguer anomaly, (D) crustal thickness, and 
(E) Bouguer gravity gradient of the Orientale basin and surroundings. Dashed lines in (A) 
from innermost to outermost correspond to the Inner Depression (ID), Inner Rook ring 
(IRR), Outer Rook ring (ORR), and Cordillera ring (CR). The solid white line in (D) shows 
the location of the cross-sectional profile A-A’ in Fig. 2a. Blue lines show the locations of 
the azimuthally averaged cross-sections in Fig. 2b. Topography is updated from Lunar 
Observer Laser Altimeter (LOLA) (77), map LDEM_64, 0.015625° spatial resolution. To 
highlight short-wavelength structure, we have subtracted spherical harmonic degrees and 
orders less than 6 from the Bouguer gravity field. Calculation of crustal thickness and 
Bouguer gravity gradient are discussed in the SOM. 


(A) Cross-section of the subsurface structure of the Orientale basin along the profile 
shown in Fig. 1D, from southeast to northwest. Crust is shown as tan, melt sheet as red, 
and mantle as green. Arrows above the cross-section denote, inward to outward, 
Orientale’s Inner Depression (ID), Inner Rook ring (IRR), Outer Rook ring (ORR), and 
Cordillera ring (CR). The heavy solid line indicates the base of the crust in the presence 


of a melt sheet that is 10 km thick, 350 km in diameter, and 2650 kg m* in density; the 
thin solid line indicates the base of the crust if the melt sheet density is identical to that of 
the crust. Vertical exaggeration is 5:1. (B) As in (A), with no vertical exaggeration and 3x 
higher resolution filter for downward continuation, azimuthally averaged in sectors 
(azimuth measured clockwise from north, see Fig. 1D) to suppress noise. Red lines give 
the location of proposed faults dipping inward at 50° dip angle from the nominal surface 
position of the CRR, ORR, and IRR. Other variations in crust-mantle boundary depth 
suggest the presence of additional faults. 
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GRAIL’s Mapping Mission 

The GRAIL twin spacecraft, Ebb and Flow, orbited in tandem in polar orbit above the 
Moon and mapped the lunar gravity field at unprecedented resolution and accuracy. 
After insertion into lunar orbit on December 31, 2011 (Ebb), and January 1, 2012 (Flow), 
the spacecraft collectively executed a succession of 19 maneuvers to achieve the science 
orbit and mapping configuration. GRAIL’s Primary Mission (PM) (39), which 
commenced on March 1, 2012, consisted of 89-days of mapping over three 27.3-d global 
mapping cycles at an orbit inclination of 89.9° and a mean altitude of 55 km. In practice 
the orbits deviated from circular due to perturbations from the lunar gravity field. When 
the solar orientation became unfavorable for satellite-to-satellite ranging on May 30, 
2012, the spacecraft discontinued mapping and the orbits were raised to a mean altitude 
of 84 km. 

During GRAIL’s Extended Mission (XM), which commenced on August 30, 2012, 
the dual spacecraft orbited above the lunar surface at a mean altitude of 23.5 km (8), less 
than half the mean mapping altitude in the PM. The low-altitude XM orbit was 
maintained by executing 2-3 weekly maneuvers (40). 

Subsequent to a lunar eclipse on November 28, 2012, the mean altitude of the orbiters 
was decreased by another factor two to 11.5 km, for additional mapping (Fig. $1). The 
principal objective of this Endgame mission phase was to perform mapping of the 
Orientale basin at the highest possible resolution. The low altitude of the spacecraft and 
variation of topography, coupled with the variability of the lunar gravity field, dictated 
that the orbital altitude coverage over the basin was non-uniform (Fig. S2), with the 
lowest altitudes achieved in the eastern part of the basin and in particular over the basin 
rings, which were mapped from 2-km altitude. GRAIL’s science mapping was completed 
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on December 14, after which a series of engineering experiments was executed prior to 
deorbit on December 17, 2012. 


GRAIL Global Gravity Solutions 

The global gravity field of a planetary body is typically represented in spherical 
harmonics, which are solutions to Laplace’s equation for the gravitational potential on a 
sphere. In practice, the resolution of gravity measured from a spacecraft will depend on 
both coverage and spacecraft attitude, and so it is common for unmodeled signal to be 
present even for high-degree and -order gravity field solutions. Published gravity models 
that incorporate global data from the GRAIL PM are to degree and order 660 (spatial 
block size = 8.3 km) (4/, 42), and published models that incorporate observations from 
both the PM and XM are to degree and order 900 (spatial block size = 6.0 km) (/0, 43). 
The analysis in this paper uses a degree- and order-1200 (spatial block size = 4.5 km) 
field (44), produced in the manner of our previously published models. Due to the 
varying spacecraft altitude that was a combined consequence of the lunar topography and 
perturbations to spacecraft orbits by the uneven subsurface mass distribution, the full 
resolution of this global model is achieved only in the areas of lowest-altitude mapping. 
Orientale contains the lowest-altitude coverage of the GRAIL mission, due to a decision 
to produce the highest-possible-resolution map of this important structure. The Orientale 
map from the degree- and order-1200 field was produced at the Jet Propulsion Laboratory 
(44). 


Local Modeling Using Short-arc Analysis 

As an alternative approach to achieving the highest-resolution map permitted by the 
data, we developed a local model of the gravitational signature of Orientale that utilized 
the residual signal remaining from a gravity model to degree and order 900. We analyzed 
GRAIL’s complete XM and endgame Ka-band range rate (KBRR) tracking data from 
240° to 300° E longitude and -55° to 5° latitude. Data were analyzed in arcs of 
approximately 18-minute duration. The short-arc analysis method (9) utilizes a spherical 
harmonic representation of the 12-epoch state parameters that quantitatively characterize 
the baseline between the satellites. First used to determine temporal gravity solutions 
(45, 46) using data from the Gravity Recovery and Climate Experiment (GRACE) 
mission (47), the method as applied to the Moon has been demonstrated with GRAIL 
observations over the lunar south pole (48). 

The gravity field may be expressed with the Stokes formulation as a gravitational 
potential W: 

W=U+T (1) 
where U is the normal potential of a reference field and T is the disturbing potential 


R 
TP) =—ff, S(P,Q)Ag(Q)do (2) 
and where the integration is taken over the surface o, S(P,Q) is the Stokes-Pizetti kernel 


(49), and the gravity anomalies (Ag) may, with a spherical approximation, be written 
Ag =-=-2- (3) 
where r is the radial distance in spherical coordinates. If T is expressed in spherical 


harmonics, then the anomalies and spherical harmonics are equivalent and anomalies can 
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be expressed in spherical harmonics by multiplying by a factor (/-/), where / is the 
spherical harmonic degree. The partial derivatives of the KBRR data with respect to the 
anomalies can be obtained from the acceleration on the satellite by differentiating a 
version of equation (2) in which the integral has been replaced by a summation over the 
separate anomalies. 

Given a GRAIL gravity model to degree and order 900 (/0) as an a priori model, 
parameters that included the baseline vector pitch, baseline rate of change of vector 
magnitude, and baseline rate of change of vector pitch were adjusted using the Goddard 
Space Flight Center GEODYN II Orbit Determination and Geodetic Parameter 
Estimation Program (50). After adjusting these baseline parameters and converging the 
short arcs, GEODYN was again used to calculate partial derivatives of the KBRR data 
points with respect to the adjusted baseline parameters and selected gravitational field 
parameters. Local anomaly values (Ag,q;) were estimated with respect to the global 
model (Ag¢rgmoooa) at the center coordinates of grid points. For our final model we 
applied a grid with mixed resolution, to account for the varying altitude above 
topography. For the area between 240° and 275° E longitude we used a resolution of 
0.1667 x 0.1667 deg’, whereas for the area between 275° and 300° E longitude, which 
included the lowest-altitude passes, we used a resolution of 0.1 x 0.1 deg”. The maximum 
spatial resolution is thus 3 to 5 km. 

We applied a neighbor-smoothing constraint (5/7, 52) to the full solution 

Agfuu = AGaaj + AGcRGms00a> (4) 
which, in addition to allowing the mapping of gravity anomalies to be extended to the 
highest resolution, had the benefit of mitigating high-degree striping evident in the 
spherical harmonic solution. Supercomputers at the NASA Center for Climate Simulation 
(NCCS) at NASA’s Goddard Space Flight Center were used to compute the solutions. 

A comparison of residuals from the global spherical harmonic model with those from 
the local model is shown in Fig. S3 and illustrates the considerable improvement of the 
latter model. The average root mean square (RMS) residual of the KBRR fits over the 
Orientale region improves from 0.987 ym s" for the global model to 0.225 wm s" for the 
local model (77% improvement). The median of the fits improves from 0.629 um s"' to 
0.160 um s’ (75% improvement). 


Local and Global Map Comparison 

The local model and degree-1200 global model are independent representations of the 
gravitational signature of Orientale produced from the same GRAIL tracking data set. 
The models are of comparable resolution and, as shown in Fig. S4, they show strikingly 
similar small-scale structure, which provides high confidence in the short-wavelength 
character of the map. 


Bouguer Anomaly Field and Crustal Structure 

In order to isolate subsurface density variations, we subtracted the gravitational 
attraction of surface topography from a spherical harmonic model to degree and order 
1080 to yield a map of Bouguer gravity anomalies. The calculation assumes an upper 
crustal density of 2550 kg m® (/8) mantle density of 3220 kg m”. In the Bouguer 
anomaly plot in Fig. 1C we subtracted spherical harmonic degrees < 6 in order to 
highlight short-wavelength structure. 
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Crustal thickness was determined in a manner analogous to an earlier study using 
GRAIL gravity and LOLA topography by our group (/8). In this work we exploit the 
considerable improvement of the XM data, with more than double the resolution of the 
data in the earlier study. Downward continuation of the Bouguer anomaly to a crust- 
mantle interface whose average radius is 1703.15 km yields an average thickness of 34 
km and satisfies global constraints. The relief on this density interface produces a 
gravitational signal that matches the Bouguer anomaly up to a degree and order 
determined by a constraint on its amplitude needed to regularize the solution in Fig. 2A. 
This constraint is equivalent to a low-pass filter whose transfer function is 0.5 at degree 
80, equivalent to a spatial feature size of ~68 km. However, the signal in the anomaly 
data over Orientale exceeds formal noise estimates at all degrees up to and beyond degree 
600, and thus warrants a higher resolution in this region, at the risk of unduly amplifying 
unmodeled density variations in the near-surface crust. We explored models with 2—3 
times higher resolution, whose amplitudes approached the limiting noise level in the data. 
Spatial averaging over azimuthal sectors, as in Figure 2B, preserves the radial resolution 
of mantle relief coherent with ring structures, while suppressing signals from the 
shallower density variations in the crust. The azimuthal sectors shown in Fig. 1D are 
chosen to avoid large E-W gradients in elevation and crustal thickness. 

The range of slopes of the region of uplifted mantle that results from a crustal 
inversion depends roughly linearly on the density contrast assumed between the crust and 
mantle. However, the range depends even more on the filter used to stabilize the 
inversion, with higher slopes generally resulting when higher spherical harmonic degrees 
are passed by the filter than in the earlier inversion (/8). Features that are revealed at 
higher resolution are suggestive of and consistent with fault displacements in the 
accompanying hydrocode calculation (/4), although faults are virtually indistinguishable 
from more continuous relief and thus not necessarily required by the data. 

Fig. S5 show radial profiles of surface topography and crust-mantle boundary relief 
that highlights the azimuthal variability of surface and subsurface structure. Within the 
inner depression, azimuthal variations in thickness are +1 km. Outside the Inner Rook 
ring, azimuthal variations in crustal thickness are +5—10 km, due in large part to 
azimuthal variability in the topography. 


Effect of a High-density Melt Sheet on the Crustal Structure Model 

Our crustal structure inversion included a central sheet with higher density than the 
surrounding crust. This sheet could be representative of crustal impact melt that formed 
initially with zero porosity, an impact melt sheet that contained both mantle and crustal 
components, or a more mafic crustal material formed by differentiation of the central 
impact melt pool. 

An end-member model with no melt sheet contribution, a crustal density of 2550 kg 
m’’, and a mantle density of 3220 kg m® (/8) is consistent with some impact models in 
which cool crustal material migrates to the basin during basin collapse and covers 
initially exposed molten mantle material (/4, 53). However, because the pre-impact lunar 
crust is largely to entirely removed by the impact (32), the crust that is measured by 
GRAIL at the basin center may be dominated by the solidified impact melt sheet. 

Impact breccias and impact melt rocks from the Apollo 14, 15, 16, and 17 landing 
sites provide samples of impact rocks from the Imbrium, Nectaris, and Serenitatis basin- 
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forming impacts. The measured bulk densities of 18 such rocks range from 2440 to 2830 
kg m° (54). These bulk densities depend primarily on the amount of clastic material with 
pre-existing porosity in the breccias and the amount of vesicularity within the impact 
melt. Because of the high level of shock heating in the basin center, we do not expect a 
large abundance of clasts in the central melt sheet. The Apollo samples have an average 
melt sheet bulk density of 2720 kg m° for clast-poor impact melt rocks. Vesicularity 
should decrease with increasing depth (pressure), so it is possible that a slightly higher 
average bulk density is appropriate for the melt sheet. 

An important consideration is the effect of differences in composition between the 
Orientale melt sheet and the Apollo impact melt samples. The Apollo samples have a 
mean FeO concentration of 6.8 weight % (range 3.1 to 9.2 %) and a mean TiO, 
concentration of 0.95 % (range 0.2 to 1.5%), whereas remote sensing observations 
indicate that the Orientale melt sheet (Maunder Formation) has mean FeO and TiO, 
concentrations of 4.4% and 0.6%, respectively (55). The scaling between composition 
and grain density of Huang and Wieczorek (56) suggests that the Orientale melt sheet is 
about 2.4% less dense than the Apollo samples. If the Orientale melt sheet and the Apollo 
impact melt samples have the same average porosity, this percentage reduction in density 
can also be applied to the bulk densities, implying a bulk density of 2650 kg m® for 
Orientale’s melt sheet. 

Relative to the typical feldspathic highland crust, the higher melt sheet density 
reduces the average density difference between crust and mantle and thus requires a 
thicker crust in the basin center. For a melt sheet density of 2650 kg m”°, the minimum 
crustal thickness in the basin center is 10 km; a denser melt sheet of 2720 kg m®* 
consistent with Apollo samples yields a minimum crustal thickness of 11 km. These 
minimum crustal thickness values are 2 and 3 km greater than a model in which the melt 
sheet has the same density of surrounding feldspathic crust (2550 kg m°). If the crust 
beneath the basin floor consists of a mixture of feldspathic crust (that either was not 
excavated or that flowed into the basin center from the surroundings during collapse of 
the transient cavity) and the impact melt pool, gravity data alone cannot constrain either 
the relative proportions of these materials or the thickness of the crust beneath the basin 
floor. 


Composition of Crust Excavated by the Orientale Impact 

Modeling the structural evolution of the Orientale basin requires knowing the 
lithologies of the crust beneath the basin. In our modeling we have assumed that 
anorthositic rock dominates the entire crustal column. However, lunar sample (62) and 
remote sensing (63) data suggest that in the Procellarum KREEP Terrain [PKT, where 
KREEP is an acronym for lunar rocks enriched in potassium (K), rare earth elements 
(REE), and phosphorus (P)], of the lunar nearside, the lower crust is composed of rocks 
richer in olivine and pyroxene than the upper, plagioclase-rich crust. It is conceivable that 
this mafic lower crust is present everywhere on the Moon. 

We test whether mafic ejecta is present in the Orientale basin and its ejecta using 
remote sensing data of the region and the expected compositional signature of LKFM. 
Korotev (64) compiled bulk compositional data for samples of LKFM impact melt 
breccias from the Apollo 14-17 landing sites. Although compositions vary from site to 
site, the chief compositional characteristics are well defined: compared with the typical 
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feldspathic crust, LKFM is significantly enriched in FeO and incompatible trace elements 
(Th, U, rare earth elements). Mineral abundances also differ: LKFM is enriched in low- 
Ca pyroxene and depleted in plagioclase compared with the feldspathic highlands. 
Olivine is variable in abundance and might represent mantle materials incorporated into 
the LKFM impact melts (64). Bearing in mind that LKFM is a mixture of a KREEP 
component, anorthosite, and a small component of olivine-rich rock (64), we can use 
Korotev’s (64) compositional compilation and modal analyses of lunar samples (65) as 
guides to establish a chemical and mineralogical fingerprint of LKFM and use remote 
sensing data to search the Orientale basin structure for that signature. 

Remote sensing data exist for the concentrations of FeO, Th, and major minerals. For 
FeO concentrations we use maps of Clementine surface reflectance data (binned to 1 km 
spatial resolution) and an iron-reflectance algorithm (66). For Th we use the Lunar 
Prospector Gamma-Ray Spectrometer (GRS) dataset (67) gridded at 0.5°, about 15 km 
spatial resolution. We use mineral data (68), which are derived from radiative-transfer- 
based mineral maps (69), validated and improved by comparison with Lunar Prospector 
GRS elemental concentrations gridded to 2° bins. The minerals of interest for detecting 
mafic impact melt are plagioclase, low-Ca pyroxene, and olivine. The LKFM component 
contains, as estimated from data compilations (64)-(65), 7-11 wt% FeO, 5-10 ppm Th, 
45-50 wt% plagioclase, 25-30 wt% low-Ca pyroxene, and 5—10 wt% olivine. These 
abundances are similar to the bulk composition of the impact melt sheet hypothesized to 
be present within the Orientale basin (60). We examined this set of remote chemical and 
mineralogical data to see if the hypothesized mafic, Th-rich crustal component is present 
in the Orientale region. 

Results are shown in Fig. S6. Fig. S6A shows the topography of Orientale and its 
surroundings. The Inner Rook ring is outlined for reference and appears on all images in 
Fig. S6. The mare basalt fill in the middle of Orientale is clearly identifiable by its high 
FeO concentration (Fig. S6B). The region with lower FeO (~4 wt%) between the maria 
and the Inner Rook ring composes the Maunder Formation, thought to be impact melt 
(e.g., UZ, 5, 70). A careful study of the composition of the basin interior (55) indicates that 
the Maunder Formation has an average FeO concentration of 4.4 + 2.0 wt%. The 
Maunder Formation is low (~1.3 ppm) in Th (Fig. S6C), as are all basin deposits, except 
for the mare basalts. Thus, both FeO and Th concentrations in the exposed Orientale 
impact melt sheet (the Maunder Formation) are substantially lower than the ranges 
observed in LKFM rocks or in regions associated with the PKT. We suggest that the low 
FeO and Th concentrations rule out a significant contribution of mafic lower crustal 
materials to Orientale surface deposits. 

Taking a broad view of the composition of Orientale deposits, the impression is that 
it is composed of plagioclase-rich rock. All areas (except for the maria) contain at least 
70 wt% plagioclase (Fig. S6-D), with many areas containing more than 85%. This result 
is consistent with the abundance of pure anorthosite in numerous small areas (7/, 72). 
Olivine and orthopyroxene abundances are generally less than 10 wt% and quite variable. 
Thus, the Orientale basin deposits all appear to be typical lunar feldspathic highlands, 
with no strong evidence for ejected lower crustal materials. We conclude that the 
assumption that the crust has the density of porous anorthite (2550 kg m”*) (/8), is 
appropriate for mapping the crustal structure outside of the Inner Rook ring. Within the 
Inner Depression, where the melt sheet is expected to be thickest, the melt sheet is largely 
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covered by a thin veneer of mare basalt, and thus remote sensing does not constrain its 
composition. We therefore used Apollo sample data to estimate the density of the central 
melt sheet, as described above. 


Excavated Volume 

In determining the volume of material excavated by the Orientale impact we assumed 
a mean global crustal thickness of 34 km (/8), a crustal density of 2550 kg m® (which for 
a 2900 kg m® grain density yields 12% porosity), and a mantle density of 3220 kg m”; 
these figures collectively satisfy the observations and a 30-km crustal thickness constraint 
at the sites of the Apollo 12 and 14 seismic stations (75). 

In this model the median thickness is 42.7 km in the region surrounding Orientale 
(Fig. 1D, Fig. 2), and the base of the mantle plug extends to a depth of 43.5 km below the 
1738-km-radius topographic datum. 

An alternative model with a mean global crustal thickness of 43 km that corresponds 
to an Apollo seismic constraint of 38 km (76) results in thicker crust at Orientale and a 
lower mantle density, which would effectively increase the volume of the mantle plug by 
approximately 10%. 


Polar Projections of Orientale 

To highlight surface and subsurface deviations from axisymmetry of Orientale, we 
plot the topography and gravity field in polar coordinates about the basin center. A re- 
projection of the datasets in Fig. 1 is plotted in Fig. S8. To generate these “unwrapped” 
maps, we linearly interpolated the topography and gravity datasets on a grid of azimuth 
(measured clockwise from north) and radial distance with respect to the center of 
Orientale. The basin center (-19.3°N, 266.0°E) was chosen such that the Bouguer 
anomaly high at the center of Orientale becomes as close as possible to axisymmetric in 
this projection. Linear interpolation of the original topography and gravity datasets is 
necessary in order to smooth the data toward the center of the basin, where latitude and 
longitude points are relatively widely spaced in an azimuth and radius projection. After 
this transform, radial features become vertical, parallel lines. Axisymmetric features (e.g., 
the Cordillera and Inner and Outer Rook rings) become horizontal, parallel lines. 
Elliptical features and axisymmetric features that are misaligned with the center of 
Orientale become sinusoids. 

As with any two-dimensional Cartesian plot of circular or spherical data, these polar 
maps can distort relative sizes, shapes, and angles. As we are interested in identifying 
radial features within Orientale, we used a modified Mercator projection: 


0 
ye In (tan (90° = 5) ford >100km 


(—3.2976 x 10°*)6 +. 4.5479 ford <100km 
where @ is the angular distance from the center of Orientale (which is related to distance 
d from the center by 0 = d/ (30.334 km/degree)). At radial distances beyond 100 km 
from the center of Orientale, this projection is conformal (angles are preserved around all 
locations, and the horizontal and vertical scales around all positions are equal). Thus, the 
shapes of small features are preserved (e.g., craters are round), although the linear scale 
changes with distance from the center of Orientale. Like the polar regions of standard 
Mercator maps, this projection becomes unusable near the basin center. Thus, we use a 
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linear relationship between y and @ within 100 km, constructed such that y(@) is smooth 
and continuous at 100 km. Although this relation prevents the singularity at 9=0° found 
in the Mercator projection, the portion of the map within 100 km of the basin center is not 
conformal. 


Bouguer Gravity Gradients 

The gravity gradient quantifies the curvature of the Bouguer gravity field, which can 
be useful for identifying otherwise hidden density structures in the lunar interior (77, 78). 
The Bouguer gravity gradient is determined from the Bouguer gravity gradient tensor, 
which is a nine-component, symmetric tensor with five independent components 
constructed by taking gradients of the Bouguer gravity field: 
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where gx, Jy, , and g, are the components of the Bouguer gravity field in a local 
Cartesian reference fame. As we are primarily interested in the horizontal gravity 
gradients (Ty, Ty, Ty), we evaluate the eigenvalues of the horizontal gradient tensor, 
I,, and I,,, which represent the maximum and minimum horizontal curvature of the 
Bouguer gravity field. We follow Andrews-Hanna et al. (77) and evaluate the maximum 


amplitude horizontal gravity gradient: 
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This maximum amplitude eigenvalue (referred to as the “Bouguer gravity gradient” or 
“gravity gradient” in this paper) effectively quantifies the gravity gradient orthogonal to 
significant density structures, as long as the structure dominates the local gravity 
gradients. For this latter criterion to be true, it is necessary that the gravity gradients be 
evaluated from filtered Bouguer gravity fields to avoid “striping” along the spacecraft 
trajectory and other sources of noise. The Bouguer gravity was filtered with a low-pass 
filter applied at degree 600, corresponding to a block size of 9 km. To further smooth the 
gravity gradients, the field was calculated on a spherical reference surface at a radius of 
1740 km, with the upward continuation of the gravity anomalies to this surface having a 
smoothing effect on shorter wavelengths. At low degree and order, topography in the 
mantle (such as the mantle uplift beneath Orientale) dominates the curvature in the 
Bouguer gravity field, and thus the gravity gradient. To highlight shorter wavelengths, 
and smaller density anomalies, we filtered the Bouguer gravity field with a high-pass 
filter, removing long-wavelength (1 S 50) anomalies. It is important to note that Bouguer 
gravity gradients can produce spurious features in gravity gradient maps because of the 
highly non-linear action of evaluating gravity gradients and determining eigenvalues. We 
have avoided use of the minimum eigenvalue, which is particularly susceptible to this 
issue. 


22 


FIGURE CAPTIONS AND FIGURES 


Fig.S1. Apoapsis (green, violet) and periapsis (blue, red) altitudes of GRAIL-A (Ebb) and GRAIL- 
B (Flow) during the GRAIL mission Endgame phase. 


Fig. S2. Minimum altitude above topography of the GRAIL dual spacecraft over the Orientale 
basin and surroundings. 


Fig. S3. (A) KBRR residuals with respect to a GRAIL degree-and-order-900 gravity model (70) 
over the Orientale basin and (B) residuals after local modeling for a smoothing factor of 10° 
(48). 


Fig. S4. Comparison of free-air anomaly fields for the Orientale basin from the GRAIL mission 
using two independent approaches to achieve high resolution. (A) Local model with neighbor 
smoothing for a scaling factor of 10° as described by Goossens et al. (48), and (B) degree- 
1200 spherical harmonic solution plotted to degree and order 900. 


Fig. S5. Radial profiles of surface topography and mantle relief (in black) across the Orientale 
basin. Basin center is at far left. Radial averages and standard deviations are plotted in red. 
The average positions of the basin rings are marked with vertical dashed lines (from left to 
right): Inner Depression, Inner Rook ring, Outer Rook ring, and Cordillera ring. Vertical 
exaggeration = 5:1. 


Fig. S6. Maps of (A) topography, (B) FeO concentration, (C) Th concentration, (D) plagioclase 
abundance, (E) orthopyroxene abundance, and (F) olivine abundance for the Orientale basin 
and nearby areas. The position of the Inner Rook ring is shown in white for reference. Data 
sources: Topography updated from Lunar Observer Laser Altimeter (LOLA), map LDEM_64, 
0.015625° spatial resolution (79); FeO from Clementine 950 nm/750 nm, 1-km data converted 
to FeO with the algorithm developed by Lucey et al. (66); Th from Lunar Prospector 0.5° 
data.; mineral maps from Crites and Lucey (68). 


Fig. S7. Lunar Reconnaissance Orbiter Wide Angle Camera (LRO/WAC) image 
[WAC_GLOBAL_E000N1800_032P.img] of the Orientale basin and surroundings showing 
distribution of highland anorthosite crust (light) and maria (dark; 21). The dark ring in the 
southern part of the basin centered on the Outer Rook ring is the product of a basaltic 
pyroclastic eruption from a central vent (80). 


Fig. S8. Polar projections of (A) topography, (B) free-air gravity anomaly, (C) Bouguer gravity 
anomaly, (D) crustal thickness, and (E) Bouguer gravity gradient with respect to the center of 
the Orientale basin, encompassing longitudes 240° to 295° E and latitudes -42° to 10°N. This 
polar projection highlights annular structures, which appear horizontal, and radial structures, 
which appear as vertical fabric. The center of the basin (-19.3°N, 266.0°E) is stretched along 
the abscissa at zero radial distance, and the radial distance of 700-km corresponds to the 
distal extent of the basin. The approximate locations of the Inner Depression (radial distance 
r=170.5 km), Inner Rook ring (r240.5 km), Outer Rook ring (r=319.5 km), and Cordillera 
(r=468.5 km) are indicated by dashed lines. 
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